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Abstract 

Four of the six parameters denning the two-loop tttt scattering amplitude have been de- 
termined using Roy dispersion relations. Combining this information with the Standard 
expressions, we obtain the threshold parameters, low-energy phases and the 0(p A ) constants 
l{, l r 2 . The result l r 2 {M p ) = (1.6 ± 0.4 ± 0.9) x 10~ 3 (l 2 = 4.17 ± 0.19 ± 0.43) reproduces the 
correct D-waves but it is incompatible with existing Standard analyses of form factors 
beyond one loop. 
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1. During the last few years there has been a noticeable revival of interest in the high preci- 
sion analysis of low-energy tttt scattering flf- ||13|| . There are at least two reasons for this. First, 
it has been shown |L3], [I], |3|] and repeatedly emphasized |TJj] that the nn scattering amplitude 



in the threshold region is particularly sensitive to the strength of quark anti-quark pair con- 
densation in the QCD vacuum: the smaller the condensate, the stronger the isoscalar S"-wave 
7T7T interaction. The accurate measurement of S'-wave scattering lengths would, indeed, provide 
the first experimental evidence in favour of, or against, the standardly admitted hypothesis 
according to which the mechanism of spontaneous chiral symmetry breaking is dominated by 
the formation of a large < qq > condensate. Within QCD, this hypothesis is by no means 
a logical necessity and its experimental test might well become an important step towards a 
non-perturbative understanding of the quark-gluon dynamics. The second reason which makes 
detailed inr studies topical, is that there are two new high precision experiments currently un- 
der preparation: i) The phase shift difference Sq(E) — 5\{E) at low energies (E < 400MeV) 
will be extracted from a new X^-decay experiment |15|] performed with the KLOE detector at 



the Frascati 0-factory DA$NE (T§. ii) At CERN, the project DIRAC |T7j aims at the mea- 
surement of the lifetime of 7t + tt~ atoms to 10%, implying the determination of the combination 
of scattering lengths | a[J — <2q | with a 5% accuracy. On the theoretical side an even better 
precision can be reached by a systematic use of chiral perturbation theory [J8[ [19[ (%PT). 
The low-energy expansion of the tttt scattering amplitude A(s\t,u) starts at order 0(p 2 ) given 
by Weinberg more than 30 years ago [pOfl . Subsequently, the one-loop 0(p 4 ) contribution to 



A(s\t,u) has been calculated by Gasser and Leutwyler |2T], [19[ . It is given by four low-energy 
constants h, I2, h, h besides the (charged) pion mass M n and the decay constant F n . The 
present state of the art involves the two-loop 0(p 6 ) order and the present letter concerns this 
degree of accuracy. 



2. The 0(p 6 ) amplitude A(s\t,u) has been first given in Ref. IJ in the form 



A(s\t, u) = A KMS F(s\t, u; a, (3\ Ai, A 2 , A 3 , A 4 ) + O 



A J 'Uu 



(1) 



The function A^msf, which depends on the Mandelstam variables s, t, u and on the six pa- 
rameters a, P, Ai, . . . A 4 , is explicitly displayed in Here, p denotes the characteristic pion 
momentum and A# is the mass scale of bound states not protected by the chiral symmetry, 
Ah ~ 4:7i F n ~ 1 GeV. The result (1) holds independently of the strength of the quark con- 
densate. The latter merely shows up in the size of the constant a: for standard, large valuesS 
< qq >~ —(250 MeV) 3 one has a ~ 1 and its value increases up to a ~ 4 for |< qq >| decreas- 
ing down to zero. The parameter (3 is less sensitive to the value of the condensate, remaining 



x < qq > denotes the single flavour condensate in the SU(2) x SU(2) chiral limit m u = = at the QCD 
scale v = 1 GeV. 
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always close to unity. It has been shown |3j that the remaining four constants Ai, . . . A4 can be 



rather accurately determined from the existing tttt scattering data |22| in the intermediate en- 
ergy range 0.5 GeV < E < 1.9 GeV using the Roy dispersion relations P5| . The latter explicitly 
incorporate crossing symmetry and consequently they strongly constrain the tttt amplitude at 
low energies. Equating the perturbative formula (1) with the Roy dispersive representation in 
a whole low-energy region of the Mandelstam plane, one infers the values of Ai, . . . A4, whereas 
the paramenters a and (3 remain essentially undetermined. The resulting Aj's are almost in- 
dependent of a and (3. Here we quote and use the central values corresponding to a = 1.04, 
(3 = 1.08 §, 

Ai = (-5.7 ± 2.2) x 10- 3 , A 2 = (9.3 ± 0.5) x 10~ 3 , 
A 3 = (2.2 ±0.6) x 10" 4 , A 4 = (-1.5 ±0.12) x 10~ 4 . (2) 

The quoted errors include experimental uncertainties on tttt phase-shifts and inelasticities 
in the medium energy region and an estimate of the systematic error arising from neglected 
higher orders in the low-energy representation (1). The errors due to the uncertainty in the 
high-energy behaviour of the tttt scattering amplitude are negligible. 

3. With the constants A« determined, Eq. (1) allows one to convert new high-precision 
experimental information on low-energy tttt phase shifts and/or threshold parameters into a 
measurement of a and (3 and finally, into an experimental determination of the quantity (m u ± 
m d) < qq > (the detailed relation between a and (3 and the condensate can be found in 
Ref. [l|]). Conversely, Eq. (1) can be used to predict, for each value of the condensate, all low- 
energy observables. It is of particular importance to assess with as much accuracy as possible 
the prediction concerning the standard alternative of a large < qq > condensate. The strength 
of the < qq > condensate is conveniently described by the deviation from the Gell-Mann- 
Oakes-Renner relation, i.e. by the parameter 

m FlMl 



m 2m \< qq >\ 



-1. (3) 



Here, m = \{m u + ma) is the running quark mass and mo is a mass scale characteristic of qq 
condensation. The standard alternative of a large condensate corresponds to mo ^ A#. In this 
special case the ratio (3) can be treated as an expansion parameter, m/mo = 0(p 2 / h 2 H ) and 
the general low-energy expansion becomes the standard chiral perturbation theory (S%PT) ||19|| . 
The complete S;\PT two-loop calculation of the 7T7r-scattering amplitude has been recently com- 
pleted by Bijnens et al. @. Not surprisingly, this calculation recovers the formula (1) giving, 
in addition, the expressions of the six parameters a, (3, A 1; . . . , A 4 in terms of i) M n , F n , ii) four 
0(p 4 ) constants li(fj), (/•*)> ^s(a*) an d ^(/-O an d finally Hi) six 0{p & ) constants r[(/i), . . . , /g(/z) 
which appear in the effective lagrangian and are renormalized at a scale /i. These expressions 
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read 



a 



1 1 \ M 2 



1 + 

— 24Z5 54 + 20/5 Z£ + 5Zf + 



6239 



14 3 
-8 h - — k 2 - 13 k 3 - - k A 



19 



769 1 

331776 7T 4 ' vr 2 V 3456 576 6 1 



1 

9 



27 



+ n l 2 ~ ~ Tq % - - 11 ) + 3r[ + Ar r 2 + 4r£ - 4r 4 



Ml 
Ft 



(4) 



1+I-2L + 2ZJ-— J 

8 7T 2 



M 2 

— + 



| fei — | As 3 — 3 A; 4 — 4Z5 + 5 ZJ 5 



8911 



+ 



331776 7T 4 tt 2 
+r£ + 4r^ - Ar r A + 12r r 5 - 4r 6 r 



512 864 18 1 8 2 



9 



(5) 



Ai 



3 1 



1 



+ 



7 11 79 

- An - - k 2 - - h + 8 l\Z + j 

6 2 2 3 1 4 9216 7T 4 



3456 7T 2 
+r r 3 -r r 4 + 6r r 5 - 2r r 6 



36 vr 2 

1 + 2272 L - 2496 1\ - 2160 V 2 - 384 1\ 

Ml 



F 2 



(6) 



— L + E 
3 2 



+ 27648 vr 5 



+ 



1 , 4, 1, ^ ;r /r 1223 

- ki k 2 k 4 + Ain r d -\ - 

3 3 3 2 4 331776 vr 4 



288 vr 2 

17 + 752 L + 3840 1\ + 1536 Z£ - 1920 1\ ) + 2r^ 



M 2 

7T 



(7) 



23 fc 
m 

18 



37 
36 



fc 2 + 



277 



+ 



1 



-9600/; ]+rZ--rl 



1990656 vr 4 41472 tt 2 
1 



19 + 5368 L - 13056 1{ 



(8) 



with 



and 



A I 



5 , 25 , 3311 

— ki~\ k 2 H 

36 72 3981312 tt 4 

J r 
3 r 6 , 



+ 



10368 tt 5 



-2 - 257 L + 336 + 840 Z£ 



cZ /x 



7< 



16tt 2 



Ti 



-72 



2 1 

3 ,73 = ~2 .74 = 2 



1 1q M 2 

167T 2 /i 2 



fc j (^) = (4Zr( M )- 7i L)L 



(These expressions are obtained from the expansions of the parameters b 



(9) 
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(11) 

, b$ originally 

given in 0, which are in one-to-one correspondence with a, (3, X\, . . . , A4. We prefer to work 
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with the latter set of parameters for the reader's convenience: explicit formulae for low-energy 
observables in terms of a, (3, A, are given in Ref. whereas similar expressions in terms 
of the fcj's are at present not available in the literature). A few points are worth recalling. 
i) The parameters a, (3, Aj are //-independent. This fact, together with Eq.(lO) fixes the 
scale dependence of the low-energy constants r[(/i). ii) Eqs. (4)- (9) fix the expansion of the 
parameters a, (3, Aj in powers of M% and logM^ (and/or in powers of the quark mass m), since 
li(fj), ■ ■ ■ JKfi) and the r^(fi) are quark mass independent. Contributions of successive chiral 
orders to a, (3, Aj can be identified by counting the powers of M%/F%. Notice that a and j3 start 
by an order 0(p 2 ) contribution (a — 1,(3 — 1) followed by 0(p 4 ) and 0(p 6 ) corrections. The 
expansions of Ai,A2 consist of 0(p 4 ) and 0(p 6 ) contributions, whereas A3 and A4 are entirely 
of order 0(p 6 ). Hi) The 0(p 4 ) constants I3 and 1 4 belong to the explicit symmetry breaking 
sector of the effective lagrangian. They represent the fine tuning of the < qq > condensate to 
its presumed large value: in S^PT, the deviation from the Gell-Mann-Oakes-Renner relation 
(3) is given by |19| 



rn 
m 



2l r 3 (fi) + 2l r M-^L 



2 



+ .... (12) 



Similarly, l\ controls the deviation of (3 from 1. On the other hand, the Aj's are independent 
of Z3 ( and only very weakly dependent on Z£) reflecting the fact that they are only marginally 
sensitive to the size of the < qq > condensate. In the sequel, we complete our definition of the 
standard X-PT by adopting the standardly used central values of l r 3 and l\ []B5| ^] : 

1\{M P ) = 0.82 x 1(T 3 , ll(M p ) = 5.6 x 1(T 3 . (13) 

Finally, the constants l\ and V 2 do not describe explicit symmetry breaking effects (they 
are coefficient of four-derivative terms in C^) and they are insensitive to the size of the quark 
condensate. They control the parameters Ai and A2. 

4. Equations (4)-(9) can be used to predict the parameters a, j3, Ai, . . . , A 4 and conse- 
quently, all low-energy 7T7t scattering observables, provided the low-energy constants h, . . . ,h 
and 7% . . . , rg are determined from the analysis of different processes. This is a path advocated 
by the authors of Ref. p|. In the present letter this kind of analysis will be confronted with 
additional expenmerataHnformation contained in Eq. (2). Bijnens et al. have used the val- 
ues (13) for l 3 and Z4; for l[ and V 2 they have taken the central values obtained from the S%PT 



analysis of Kn form factors ||25|| : 



l[(M p ) = -5.40 x 10-\ l r 2 {M p ) = 5.67 x 10~ d . (14) 

As for the 0{p 6 ) constants ?"[(//), the authors of take r[(lGeV) = and they check that 
this approximation confronted with a resonance saturation model produces a negligible error. 
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With the values (14), and r\ — at [A — 1 GeV one obtains (in this letter we always use 
F n = 93.2 MeV and M T = 139.6 MeV): 



a = 1.074 , p = 1.105 , 
Ai = -8.91 x 10~ 3 , A 2 = 14.5 x 10~ 3 , (15) 
A 3 = 2.04 x 10" 4 , A 4 = -1.79 x 10~ 4 . 



20 




E fMeV) 



Figure 1: The phase shift difference 5® — 5{ in the energy region of decays. The dashed 
curve is obtained with the values of Eq. (15) and it coincides with the curve displayed in Fig. 1 
of Ref. H/. The solid line is obtained with the values of Eqs. (2) and (18) while the shaded area 
results adding the corresponding error bars quadratically. The experimental points are from 
Ref. 

For these values of the parameters a, [3, X\, . . . , A4, one obtains the 5*-wave scattering lengths 
aj] = 0.218, aj] — <2q = 0.259 corresponding! to the predictions given in Ref. 0. The resulting 
phase shift difference <5q — S\ (measurable in decays) is shown as a function of the center of 
mass energy as the dashed line in Fig. 1. It reproduces the curve displayed in Fig. 1 of Ref. 0. 
Finally, a few remaining threshold parameters not discussed in Ref. are collected in the first 
column of Table 1, using the expressions displayed in Appendix D of Ref. 0]. 

2 Actually these have to be compared with the numbers given in Eq. (4) of Ref. in parentheses GeV) = 
0). The small difference provides an estimate of 0(p 8 ) effects: it is entirely due to the fact that the amplitude 
Akmsf of Ref. jl| coincides with the amplitude calculated in Ref. only modulo 0(p s ) contributions. 
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Bijnens et al. [[| 


KMSF 


Experiment [26] 


"o 


0.218 


(0.2156) 


0.209 ± 0.004 


0.26 ±0.05 


b° 


0.273 


(0.271) 


0.255 ±0.010 


0.25 ±0.03 


-10ag 


0.411 


(0.4094) 


0.44 ±0.01 


0.28 ±0.12 


-106g 


0.709 


(0.704) 


0.80 ±0.02 


0.82 ±0.08 


a o ~ a o 


0.259 


(0.2565) 


0.254 ±0.004 


0.29 ±0.05 


lOaJ 


0.395 


(0.3956) 


0.373 ±0.008 


0.38 ±0.02 


I0 2 b{ 


0.785 


(0.784) 


0.55 ±0.07 




10 2 a° 


0.263 


(0.267) 


0.16 ±0.02 


0.17 ±0.03 


10 3 a| 


0.237 


(0.2356) 


0.09 ±0.13 


0.13 ±0.30 


10 4 a 4 


0.428 


(0.478) 


0.49 ±0.07 


0.6 ±0.2 



Table 1: Threshold parameters of tttt scattering (in units of M v +) in the standard framework 
using the two-loop expressions of Ref. /[J/, App. D. The first column results from the values of 
Eq. (15) (see the text for the numbers in parentheses). The second column is obtained in the 
same way but taking the values of Eqs. (2) and (18) as input. 



The numbers in parentheses are obtained keeping in higher orders only those components 
of a, (3, Ai and A2 that actually contribute at most to the order 0(p 6 ). These exactly coincide 
with the corresponding predictions one would obtain using the amplitude given in ||. Among 
the latter it is worth noticing the value predicted for the isoscalar D-wave scattering length 
a° = 26.3 x 10~ 4 , which is three standard deviations above the value extracted from the analysis 
of Roy equations [^6fl . This disagreement reflects the fact that the value (15) of A2 is significantly 
above the value (2) inferred from experimental phase shifts in Ref. ||. We would like to stress 
that both the canonical value a° = (17 ± 3) x 10~ 4 and the determination of the constant 
A2 = (9.3 ±0.5) x 10~ 3 are based on the Roy dispersion relations |2l| using the experimental nn 
data above 500 MeV as input. Furthermore, in both cases, the dominant contribution comes 
from the P-wave in the p(770) region, which explains the relatively small error bars. These 
facts suggest that the predictions of Ref. based on (15) have to be revised in order to agree 
with the values (2) of the parameters Ai, . . . , A 4 and with the standard value of a®. We therefore 
proceed as follows: fixing l r z and l\ according to Eq. (13), we solve Eqs. (6) and (7) for /[(M p ), 
W), 



ll(M p ) = (-4.0 ± 1.0) x 10~ 3 
F 2 (M p ) = (1.6 ±0.4) x 10" 3 ± 



•l.lr^±1.0r^-6.3r^±2.1r; 
0.1r5-3.5rl±0.5rr-0.2r f 



(16) 

M=lGeV 



(17) 

/i=lGeV 



where the values and errors (2) have been used for Ai and A2. Eqs. (16) and (17) are then 
inserted back into the formulae (4) and (5) for a and (3. Keeping in mind that a and f3 are 
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sensitive to l\ and I2 only at next-to-next-to-leading level, the unknown constants r[(l GeV) 
are viewed as a source of uncertainty in a and (3. Inspired by naive dimensional analysis [01 
we take in the expressions for a and /3, r[(l GeV) = (0 ± 2) x 10~ 4 . Adding the corresponding 
uncertainties quadratically, we obtain 

a = 1.07 ±0.01 p = 1.105 ±0.015. (18) 

It should be stressed that the error in Eq. (18) does not include the uncertainty in the low- 
energy constants 1% and l\. As in the case of the chiral condensate itself, the constant l\ has not 
yet been determined experimentally and for this reason it is hard to associate an error bar with 
it. The values (18) have to be viewed as corresponding to the "standard alternative" of a large 
condensate defined by the values (13) of l\ and l\. We now use the formulae given in Ref. |lj] 
to generate the predictions for threshold parameters and phase shifts that correspond to a, (3 
(18) and Ai,...,A 4 (2). Adding the errors quadratically, the resulting threshold parameters 
are summarized in the second column of Table 1. One observes that the deviations of a[j 
and a[j — Oq from their central experimental values are significantly larger than predicted in 
Ref. [Q] . Notice that now, the Z)-wave scattering lengths perfectly agree with their Roy-equation 
"experimental" values as expected from the manner the values (2) of the constants Ai, . . . , A4 
have been obtained. 




300 400 500 



Figure 2: The isospin 2, S-wave phase shifts at low energies. The different curves are 
obtained with the values (15) (dashed) and the values (2) and (18) (solid). In the latter case 
the shaded area shows the corresponding error band. The experimental points are taken from 
Ref. fl/. 
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A similar conclusion holds for the phase shift difference 5® — 5\, shown as the solid curve 
in Fig. 1 with the error band indicated by the shaded area: the curve displayed in Ref. |2[ is 
significantly higher, i.e. closer to the experimental central-value points. For illustration, the 
phase 5q is also shown in Fig. 2. 

5. We finally address the question of interpreting the mismatch described in the previous 
paragraph. Its origin clearly appears upon comparing eqs (16) and (17) with the values of the 
constants l\ 2 (M p ) extracted in Ref. from the "unitarized" one- loop S%PT Kn form factors 



(Eq. (5.10) of p5]). Including errors the latter read: 

l[(M p ) = (-5.4 ±1.1) x 10" 3 , l r 2 (M p ) = (5.7 ± 1.1) x 1(T 3 . (19) 

The question is how close the expressions (16) and (17) can be brought to these values keeping 
at the same time the 0(p 6 ) constants 7*3(1 GeV), . . . , r£(l GeV) at a reasonable size. If one 
proceeds as before treating the r['s at 1 GeV as randomly distributed around zero with a 
standard deviation ±2 x 10~ 4 , one gets: 

l[(M p ) = (-4.0 ±1.0 ±1.8) x 1(T 3 , 

l r 2 (M p ) = (1.6 ±0.4 ±0.9) x 10" 3 , (20) 

or 

J x = -0.37 ±0.95 ±1.71 , l 2 =4.17 ±0.19 ±0.43 , (21) 

where the first error has its origin in Ai and A2 (Eq. (2)), whereas the second error arises from 
the presumed uncertainties in the individual r^s added quadratically. Two cheks of the size of 
the constants r, are conceivable, i) First, one can make a full use of informations contained in 
Eq. (2) determining the parameters l[ 2 (M p ) and rj(l GeV), . . . ,rg(l GeV) by a simultaneous 
fit to Eqs. (6)-(9) and to the constraints r[(l GeV) = ± 2 x 10~ 4 . The resulting x 2 /d.o.f. is 
1.9/2 and one obtains 

ZJ(Mp) = (-4.0 ± 0.5) x 10" 3 , l r 2 (M p ) = (2.0 ± 0.3) x 10~ 3 , (22) 

compatible with (20), whereas for the r,'s one gets 

7-5(1 GeV) = (-0.3 ± 2.0) x 10" 4 , rj(l GeV) = (-0.7 ± 0.9) x 10" 4 , 



(23) 



r 5 r (l GeV) = (1.5 ± 0.5) x HT 4 , r r 6 (l GeV) = (0.4 ± 0.9) x 10" 



This result turns out to be rather stable: if one increases the uncertainties in the r^s to 
±3 x 10~ 4 , the new x 2 '/d.o.f. = 1.24/2, the values (22) become (-4.8 ± 0.5) x 10~ 3 and 
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(2.1 ± 0.3) x 10~ 3 respectively, and the changes in the r^'s also remain rather modest. On the 
other hand, the errors obtained by this procedure (increase of \ 2 by one unit) and shown in 
Eqs. (22) and (23) are probably heavily underestimated, ii) Next, it is instructive to confront 
the previous discussion with the estimate of the constants r« by resonance saturation as quoted 



recently by Hannah [27] : 



r x = -0.61 x 10" 4 , r 2 = 1.3 x 10" 4 , 

r 3 = -1.70 x 10" 4 , r 4 = -1.0 x 10" 4 , (24) 
r 5 = 1.14 x 10" 4 , r 6 = 0.3 x 10~ 4 . 

Estimating low-energy constants by resonance saturation does not, in principle, fix the renor- 
malization scale n at which the estimate is supposed to hold. Actually, if a constant exhibits a 
strong scale dependence, its resonance saturation estimate is subject to caution. Interpreting 
Eqs. (24) as values of r[(/i) at /x = 1 GeV, one observes a striking coherence with the preceding 
analysis: (24) is, indeed, consistent not only with dimensional analysis or with the assump- 
tion |r[| < 2 x 10~ 4 but, moreover it agrees with the fit (23). One can even repeat the fit to 
Eqs. (6)-(9) constraining r[(l GeV) to the values (24) allowing for a 100% error: the fit is excel- 
lent (x 2 /d.o.f. = 0.91/2) and it yields 1\{M P ) = (-4.0±0.5) x lO" 3 , l r 2 (M p ) = (2.1 ±0.3) x lO" 3 , 
again compatible with (20) and (22). On the other hand, one finds that between /j = M p and 
H = 1 GeV, only the constants r 4 , r 5 and r 6 show a moderate scale dependence: had we 
assumed that the values (24) concern the scale \i = M p (as suggested in Ref. [23), the compar- 



ison with our previous analysis would be less favourable as far as the constant r 3 is concerned, 
7*3(1 GeV) = —4.9 x 10 -4 in this case. Notice however that according to Eq. (17) the correction 
to the "critical" constant l 2 (M p ) is dominated by r\ whose scale dependence is rather weak: 

r\{\ GeV) = r\{M p ) - 7 x 10" 6 . (25) 

In order that the constant l 2 (M p ) (17) differ from the value (19) by at most two standard 
deviations, the constant 7*4(1 GeV) would have to be 7*4(1 GeV) ~ — 5 x 10 -4 . This cannot be 
excluded but it looks unlikely in the light of the present analysis. 

6. The constants l\ and l 2 (19) have not been obtained from a full two-loop analysis of 
K14 form factors F and G, which is not yet available. Instead, their determination is based on 
matching a dispersive representation for the form factor F with the one-loop S^PT expressions, 
the latter merely serving to fix the subtraction constants. This method of "improving" one-loop 
%PT calculations has been often used in the past [^8| and it suffers from a basic ambiguity: one 



has to assume that the one-loop and two-loop amplitudes practically coincide in a particular 
kinematical point M. Even if one admits the very existence of such a matching point M, the 



results can still depend on its choice. In Ref. p5| the matching point has been chosen at the 
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threshold s n = 4M% of the S- wave amplitude mi —>■ K + axial current, where s n stands for the 
dipion invariant mass squared. We have repeated the analysis of Ref. [^5] for other choices of the 



matching point between s n = 4M% and the left-hand-cut branch point s n = 0. We reproduce 
the result (19) and find that it is actually rather insensitive to the matching point except in 
the vicinity of the singular point s n = 0, where the outcome for l\ (but not l 2 ) becomes less 
stable. For instance, with the matching point at s n = 2M%, we obtain 

ZJ(Mp) = (-4.8 ± 2.1) x 1(T 3 , l r 2 {M p ) = (5.3 ± 1.0) x 10" 3 . (26) 

Given the present state and quality of K14 experimental data, it seems hard to ascribe the 



discrepancy described above to the S%PT analysis performed in Ref. p5| . On the other hand, 
it should be kept in mind that outside the standard framework, i.e. for low values of the 
condensate < qq >, the constants l\ and I2 extracted from K14 data will be modified already 
at the one- loop level: since in G%PT the loop contributions are more important, the resulting 
central values of |Zi| and | Z2 1 are expected to come out somewhat smaller f2T 



7. A few concluding remarks are in order. The past determinations |19|, [H], |6|, [7| of the 
constants l\ and l 2 have operated within the 0(p 4 ) order of ;\PT. They have shown an apparent 
coherence and compatibility with the Ku analyses of Ref. P5fl . This compatibility might be 
lost at 0(p 6 ) order and we have to understand why. The resonance saturation models are the 
only ones that determine the constants directly, integrating out the resonance degrees of 
freedom from an extended effective lagrangian C e ff. However, incorporating resonances into 
C e ff is not free of ambiguities, especially if one aims at the 0(p 6 ) accuracy. On the other 
hand, less model-dependent sources of information, such as tttt D-waves [^9j and/or sum rules 
H primarily determine the physical parameters Ax, A2. It turns out that this determination 
is rather stable and barely affected by switching from order 0(p 4 ) to 0(p 6 ). At the 0(p 4 ) 
level, i.e. neglecting in Eq. (1) the two- loop effects and setting A3 = A4 = 0, one would get 
from the a° and a| experimental central values Ai = —6.4 x 10~ 3 and A2 = 10.8 x 10 -3 , to be 
compared with Eq. (2). In other words, the relationship between D-wave scattering lengths and 
the parameters Ai, A2 is almost unaffected by 0(p 6 ) effects. The latter however become rather 
important in the relationship between A2 and l r 2 . Rewriting Eq. (7) to make the dependence on 
l\{Mpj appear explicitly, one obtains 

A 2 = {l r 2 {M p ) + 5.45 x 10~ 3 } + {0.32 x l r 2 (M p ) + 1.7 x KT 3 } (27) 

where the first (second) curly brackets collect all 0(p 4 ) (0(p 6 )) contributions (r^ has been 
neglected). The 0{p 6 ) contribution is as large as 30% and it is dominated by double logs, 
whose importance has been anticipated by Colangelo 0. It follows that for a given A2 (D- 
waves), the resulting value of l 2 (M p ) can easily differ by a factor ~ 2 depending whether in 
Eq. (27) one includes the 0(p 6 ) term or not. Whether the consistency with form factors can 
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be understood within the large condensate hypothesis remains to be clarified. It might be, for 
instance, that at 0(p 6 ) level the K14 form factors also receive an important contribution from 
double logs, which the unitarization procedure would not take into account . Independently 
of this issue, the main conclusion of this letter is the following: the predictions of S%PT for a[j, 
Oq — (Zq and Sq — 8\ given in Ref. || are systematically overestimated as shown in Fig. 1 and 
Table 1 of the present paper. A closely related fact is the failure of the values of l[ and V 2 used 
in Ref. to describe the .D-waves in agreement with Roy equations analyses. This agreement 
is nicely recovered if instead the present determinations of Eq. (20) are used. This shows, 
once more, that a sensible and sensitive test of QCD in low-energy tttt scattering should be 
based on a global analysis making use of all theoretical constraints and all pertinent low-energy 
observables. 
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